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ABSTRACT 

Integrated spectroscopy is the method of choice for deriving the ages of unresolved stellar systems. 
However, hot stellar evolutionary stages, such as hot horizontal branch stars and blue straggler stars 
(BSSs), can affect the integrated ages measured using Balmer lines. Such hot, "non-canonical" stars 
may lead to overestimations of the temperature of the main sequence turn-off, and therefore underes- 
timations of the integrated age of a stellar population. Using an optimized H/3 index in conjunction 
with HST/WFPC2 color-magnitude diagrams (CMDs), we show that Galactic globular clusters ex- 
hibit a large scatter in their apparent "spectroscopic" ages, which does not correspond to that in their 
CMD-derived ages. We find for the first time that the specific frequency of BSSs, defined within the 
same aperture as the integrated spectra, shows a clear correspondence with H/3 in the sense that, at 
fixed metallicity, higher BSS ratios lead to younger apparent spectroscopic ages. Thus, the specific 
frequency of BSSs in globular clusters sets a fundamental limit on the accuracy for which spectroscopic 
ages can be determined for globular clusters, and maybe for other stellar systems like galaxies. The 
observational implications of this result are discussed. 

Subject headings: blue stragglers — globular clusters: general — galaxies: star clusters — galaxies: 
stellar content 



1. INTRODUCTION 

A common method for estimating the ages of unre- 
solved stellar systems is to measure Balmer lines and 
metal lines from integrated spectra, and compare them 
to stellar population models. The method relies on the 
fact that Balmer lines are mostly sensitive to the effective 
temperature (T e ff) of the main sequence (MS) turn-off of 
a stellar population (e.g. Buzzoni et al. 1994). This, in 
turn, provides a measure of the population age. 

However, a longstanding concern over the use of 
Balmer lines to estimate spectroscopic ages is the effect 
of non-canonical evolutionary stages on the integrated 
stellar population spectra. Particularly worrisome are 
the possible effects that bluer horizontal branch (HB) 
stars and blue straggler stars (BSSs) may have on the 
inferred T c g of the MS turn-off. HB stars are helium- 
burning stars which occupy a region in the CMD with 
a typical absolute V magnitude of ~ 0.7 mag (e.g. Har- 
ris 2001). BSSs are identified as blue, luminous exten- 
sions of MS stars (Sandage 1953). They are thought to 
form due to H refuelling processes after the MS stage, 
either from collisional stellar encounters (e.g. Davies et 
al. 1994) or from mass-transfer binaries (McCrea 1964). 
Both HB stars and BSSs with T e g > 6500 K show promi- 
nent Balmer lines in their spectra, which can mimic the 
presence of younger stellar populations (e.g. Rose 1985, 
Lee et al. 2000, Schiavon et al. 2004; Trager et al. 2005). 

In a recent paper, Cervantes & Vazdekis (2008) defined 
an optimized line-strength index for H/3, called H/3 Q , that 
minimizes the metallicity dependence of the Balmer line 
in favor of its age sensitivity. Interestingly, as pointed out 
in that paper, the integrated spectra of Galactic globu- 
lar clusters (GGCs) from Schiavon et al. (2005; here- 
after S05) exhibit a clear intrinsic scatter in their H/3 Q 
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strengths, particularly with increasing metallicity. 

In order to shed light on the above issue, in this let- 
ter we combine resolved and unresolved data of such 
GGCs to investigate whether CMD-based age differences 
among GGCs and/or different relative contributions of 
non-canonical stages can be responsible for the distinct 
integrated H/3 Q indices. In Section [5] we discuss the data 
used in this study, the analysis of which is described in 
Section [3] Finally, in Section [3] we present our conclu- 
sions and discuss our findings within the context of esti- 
mating ages for unresolved stellar populations. 

2. THE DATA 

Integrated optical spectra for 41 GGCs were taken 
from S05. These data were obtained by drift-scanning 
the core diameter of each GGC -with a spectroscopic 
aperture equal to this diameter- in order to construct a 
representative integrated spectrum. The spectra cover 
a wavelength range of ~ 3350 - 6430A with a FWHM 
of ~ 3.1A and typical signal-to-noise ratios (S/N) of 
~ 250 A^ 1 in the continuum of the H/3 line. We refer 
the reader to S05 for more details on these data. 

In Figure Q] we present the indices H/3 , H/3lick 
(Worthey et al. 1994) and [MgFe] (Gonzalez 1993) mea- 
sured for the 41 GGC spectra at the S05 spectral resolu- 
tion. Uncertainties in the index measurements (la error 
bars) account for the S/N spectra and the typical radial 
velocity error provided by S05 for each GGC. To guide 
the eye, based on the MILES stellar library (Sanchez- 
Blazquez et al. 2006a; Cenarro et al. 2007a), an extension 
of the simple stellar population (SSP) models in Vazdekis 
(1999) are overplotted at 3.lA spectral resolution (here- 
after V99+). The fact that most GGCs lie below the 
model grids arises from the well-known zero-point prob- 
lem of SSP models (e.g. Vazdekis et al. 2001; Schiavon et 
al. 2002), although this does not affect our results which 
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are based on relative differences. 

As expected, in the H/3 plot (Fig. Q}i) GGCs follow - 
on average- an old sequence in mctallicity. However, by 
simple visual inspection they seem to separate into two 
groups, particularly at the high metallicity end ([MgFe] 
> 1.2). Hence, solid and open symbols are employed 
to indicate, respectively, GGCs with high and low H/3 
values at fixed [MgFe] . Henceforth, these are referred to 
as GGCff and GGCL respectively. 

When using the H/3lick definition (Fig. [lb), both 
GGC groups are still distinguished, although the greater 
age-metallicity degeneracy of H/?lick would make their 
differentiation less clear to detect if the distinct symbol 
codes were not present. Even so, Puzia et al. (2002) al- 
ready pointed out the existence of unexpected H/3lick 
differences among certain metal rich (MR) GGCs. It 
therefore appears that some property differs between 
both groups leading to different Balmer-line strengths 
at a given metallicity. In fact, since the effect seems 
to increase with the increasing metallicity, we focus our 
analysis on those 25 GGCs from S05 with [Fe/H] > -1.35 
([MgFe] > 1.2) among which clear differences in H/3 at 
fixed metallicity are observed. 

Together with the indices in Fig. [TJ the adopted CMD- 
derived parameters for the 25 GGCs are listed in Table[TJ 
The HB morphology, as measured by the HBR parameter 
(Lee et al. 1994), and the specific frequency of RR Lyrae 
variables, SRRLyrae> are taken from the Harris (1996) cat- 
alogue (Feb 2003 revision; hereafter H03). We adopted 
the relative age estimates from De Angeli et al. (2005) 
and Recio-Blanco et al. (2006), whose applied the so- 
called vertical method on the GGC HST/WFPC2 snap- 
shot catalogue of Piotto et al. (2002). Also based on that 
catalogue, Recio-Blanco et al. (2006) estimated the max- 
imum T e g along the HB, T c gHB -considered as an HB 
morphology parameter-, whereas Moretti et al. (2008; 
hereafter M08) measured the logarithm of the number of 
BSSs inside the core radius, r c , normalised to the sam- 
pled luminosity -in units of 10 4 L0- in the F555W HST 
band within the same aperture. The last quantity can be 
considered as a logarithmic specific frequency of BSSs in- 
side the GGC r c , hereafter 5g c ss , and is representative of 
the spectroscopic data in S05 as they are both computed 
within the same aperture. 

3. ANALYSIS 

With the aim of constraining the origin of the intrinsic 
scatter in the Balmer line-strengths of our GGC subsam- 
ple, in Figured] we show the CMD-derived parameters of 
Section [2] -where available- as a function of the GGC 
metallicity from H03. Symbol codes are kept as in FigQ] 
except for the MR GGGffs NGC6388 and NGC6441, 
which being well-known second parameter clusters (Rich 
et al. 1997), are plotted as open stars -rather than solid 
squares- to facilitate further discussion. 

In Fig. EK we see that there is no dependence of the 
group location on the CMD-derived age. The difference 
in H/3 strengths between both groups is therefore not due 
to age differences among the GGCs, as would normally be 
inferred from a classical SSP index-index analysis. Note 
also that the typical dispersion of ~ 1 Gyr quoted by 
De Angeli et al. (2005) among the CMD-derived ages of 
intermediate metallicity GGCs could never explain the 
large scatter in H/3. We also rule out the possibility that 



the number of RR lyraes in the instability strip is playing 
a significant role, as the two GGC groups are well mixed 
(Fig. [2b). In addition we find that there is no obvious 
dependence on the HB morphology, as measured by ei- 
ther the HBR parameter or the maximum T e ff of the HB 
(Fig. [2b and Fig. [2ji, respectively). As expected, the sec- 
ond parameter clusters NGC6388 and NGC6441 (open 
stars) do stand out of the general trends in both panels. 
Their high H/3 strengths are naturally explained by the 
addition of hot HB stars in their integrated spectra. 

Besides basic properties of the GGC CMDs, literature 
estimates of [a/Fe] for the GGCs in S05 (where avail- 
able) show a high degree of homogeneity (e.g. Pritzl et 
al. 2005), so differing levels of a-elements cannot account 
for the observed differences. We have also ruled out that 
the Balmer lines of GGCLs are sistematically filled-in 
by emission. In fact, emission was found and corrected 
by S05 for only NGC6171 and NGC6553 (GGGffs), 
and NGC6352 (GGCL). Interestingly, the fact that 
NGC6352 has the weakest H/3 line of the sample may 
suggest that a residual emission could still be present. 

Having rejected the above mechanisms from being re- 
sponsible for the observed differences in H/3 between 
GGC-ffs and GGCLs -except in the obvious case of 
NGC6388 and NGC6441-, in Fig [2> we show the GGC 
metallicities versus S^gg. Interestingly, GGGffs and 
GGCLs -which were identified spectroscopically- sepa- 
rate cleanly into two groups in terms of their BSS specific 
frequencies. At a given metallicity, GGCs with higher 
Sggg values exhibit stronger H/3 lines, suggesting that 
BSSs are indeed affecting their integrated spectra. 

To reinforce this result, we quantify the impact of BSSs 
on the integrated spectrum of NGC6342, the GGCff 
with the highest Sggg. Based on photometric data from 
H03 and M08 for this GGC and its BSS population, we 
obtain that 13% of the GGC flux in V band within r c 
comes from a population of 7 BSSs with 0.22 < B — V < 
0.52, and a luminosity-weighted B — V of 0.33. For each 
BSS, assuming [Fe/H] = —0.65 and its B — V, we es- 
timate T e ff using the (B — V) — T g relation for dwarfs 
from Alonso et al. (1996). The T cS values of the 7 BSSs 
are in the range 5925 < T cS < 7500 K. We then com- 
pute spectral BSS templates with the above parameters 
on the basis of the MILES stellar library, using the in- 
terpolating algorithm in Vazdekis et al. (2003; Appendix 
B). After the 7 BSS templates are scaled according to 
their V fluxes and subtracted from the NGC6342 spec- 
trum, its integrated H/3 and H/3lick indices decrease by 
0.70 A and 0.62 A respectively. Taking into account that 
the averaged 5g SS of GGCLs with [Fe/H]> -0.8 (simi- 
lar to that of NGC6342) is ~ 1.4 (Fig. Eh), and assuming 
similar T c ff values for their BSS populations, the relative 
offsets in H/3 Q and H/3lick between NGC6342 and the 
above GGCLs are, respectively, 0.58 A and 0.51 A. This 
agrees with the differences between GGGffs and GGCLs 
obtained at this metallicity regime ([MgFe]> 2; Fig. []}, 
supporting the idea that BSSs are responsible for the ob- 
served differences in H/3. Even more, the effect of BSSs 
is also detected among the subsample of GGCLs with 
[Fe/H]> -0.8. NGC0104 and NGC6624, with extreme 
values of Sggg, pose the lowest and largest H/3 values 
respectively. 

4. DISCUSSION 
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Based on the close correspondence between the spe- 
cific frequency of BSSs in GGCs with [Fe/H] > -1.35 
and their integrated H/3 strengths at fixed metallicity, we 
conclude that BSSs are primarily responsible for the H/3 
variations observed in the integrated spectra of GGCs of 
intermediate-to-high metallicity. Far from discussing on 
the origin for the distinct S^ss vames among GGCs (see 
M08 for a thorough study on this topic), we here ana- 
lyze the implications of the above result in the context 
of age-dating unresolved stellar populations. 

First, caution must be employed in Balmer-line based 
age-metallicity studies of unresolved extragalactic glob- 
ular clusters (EGCs). Cenarro et al. (2007b) already re- 
ported the existence of EGCs with strong Balmer lines 
that were consistent with hosting an additional popu- 
lation of either blue HB stars and/or BSSs. Since the 
BSS fraction of EGCs is generally not known, the find- 
ing in this letter sets a fundamental limit to the reli- 
ability with which ages may be determined for EGCs 
using Balmer lines and SSP models. Taking the S05 
data as a representative old GC system, we can esti- 
mate this limit from the averaged offsets in the mea- 
sured H/3 lines of GGCiJs and GGGLs. Since the off- 
sets seem to vary with metallicity, local linear fits to all 
GGCffs and GGGLs in Fig.[T]with [MgFe]> 1.5 ([Fe/H] 
> —1.0) have been performed. For instance, at the lo- 
cation of 47Tuc (NGC0104; [MgFe]~ 2.31), we obtain 
AH/3 G = 0.46 ± 0.03A and AH/? L ick = 0.33 ± 0.04A 
(at the S05 resolution) with uncertainties accounting for 
the standard errors of the means. Thus, assuming that 
GGGLs are ~ 14 Gyr old -the largest SSP age in Fig. [TJ- 
the two offsets can be consistently misinterpreted on the 
basis of SSP models as GGC-ffs being ~ 6 — 7 Gyr old, 
that is, as a rejuvenation of up to ~ 8 Gyr. Differences 
between GGCffs and GGCLs also exist for the Lick H7 
and H<5 indices (Worthey & Ottaviani 1997), although 
they are not so apparent, probably due to their limited 
age-disentangling power for old SSPs. For these indices, 
the above test leads to rejuvenations of up to ~ 4 — 5 Gyr. 

The role of metallicity in the present discussion is 
worthwhile considering. Although the relation between 
Sggg and metallicity is not statistically significant - 
but marginally positive- over the entire GGC sample 
(in agreement with M08), we find clear correlations for 
GGCffs and GGGLs separately, as illustrated by the 
solid lines in Fig. [2^. The different slopes seem to indi- 
cate that, when BSSs are important, their relative con- 
tributions are larger at high metallicities. Note that the 
fading with metallicity expected in F555W band for the 
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most MR GGCs only accounts for up to ASggg ~ 0.2, 
so the above trends are irrespective of this effect. 

These results may also have important consequences 
for EGC studies. To understand the origin of color 
bimodality in GC systems within a context of galaxy 
formation, age-dating GC subpopulations -through the 
analysis of their integrated Balmer lines- is a com- 
mon practice (see Brodie & Strader 2006 and references 
therein). Interestingly, some papers have reported that 
the MR GC subpopulation of certain galaxies show on 
average a smaller mean age and a larger age scatter 
than their metal poor (MP) counterparts (e.g. Puzia et 
al. 2005). Although the present finding does not rule out 
the existence of true age differences between MP and MR 
GC subpopulations, the increasing importance of BSSs 
with metallicity might, at least, partially affect the re- 
sults of previous work. 

Whether all the above results can compromise the in- 
tegrated ages of other stellar systems, like galaxies, may 
rely on the mechanism that dominates the formation of 
BSSs. If stellar encounters were driving the BSS popula- 
tion, then one should not expect a major effect in galaxies 
because of their much lower stellar densities. However, 
this would not apply if mass-transfer binaries were the 
progenitors of most BSSs. In fact, Momany et al. (2007) 
and Mapelli et al. (2007) support the last scenario to 
explain the large BSS populations of dwarf spheroidal 
galaxies, and Han et al. (2007) have demonstrated the 
importance of binary interactions to understand the UV- 
upturn of elliptical galaxies (Es). It therefore seems that 
BSSs could play a non-negligible role in the integrated 
spectra of galaxies as long as they host an important 
fraction of binary stars. If this were the case and the 
potential increasing importance of BSSs with metallicity 
would still hold for massive Es, then BSSs could con- 
tribute to the age scatter reported for massive Es and to 
the fact that younger Es have higher metallicities than 
older Es (e.g. Trager et al. 2000; Sanchez-Blazquez et 
al. 2006b). This picture, however, requires of further in- 
vestigation which is out of the scope of this paper. 
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TABLE 1 

Galactic globular clusters in Schiavon et al. (2005) with [Fe/H] > —1.35 



GGC 


Group a 


H/3 Q b 


H/3t TCK b 


[MgFel b 


[Fe/HT 


1 IN t^ri.lVl 


HD D T , r *^ 


HBR C 


loo(Tnff HB) C 

° v y\ 611 J-A-Ljy 


qr c { 
°BSS 


NGC0104 


L 


2.30 ±0.02 


1.55 ±0.01 


2.31 ±0.02 


-0.76 


0.99 


0.2 


-0.99 


3.756 


1.03 


NGC1851 


11 


2.99 ± 0.03 


2.27 ±0.02 


1.32 ±0.02 


-1.22 


0.81 


13.5 


-0.36 


4.097 


1.48 


NGC2808 


L 


2.69 ±0.01 


2.01 ±0.01 


1.32 ±0.02 


-1.15 


0.76 


0.3 


-0.49 


4.568 


0.94 


NGC5904 


H 


3.06 ±0.01 


2.42 ± 0.01 


1.18 ±0.01 


-1.27 


0.81 


37.7 


±0.31 


4.176 


1.07 


NGC5927 


L 


2.43 ±0.03 


1.38 ±0.03 


2.87 ± 0.02 


-0.37 


0.92 


0.0 


-1.00 


3.724 


1.40 


NGC6121 


H 


3.12 ±0.03 


2.28 ± 0.03 


1.34 ±0.02 


-1.20 


0.91 


52.7 


-0.06 


- 


- 


NGC6171 


H 


2.81 ±0.05 


2.08 ± 0.05 


1.76 ±0.03 


-1.04 


0.99 


31.0 


-0.73 


3.875 


1.68 


NGC6266 


H 


2.88 ±0.02 


2.09 ± 0.02 


1.41 ±0.02 


-1.29 


0.92 


15.6 


±0.32 


4.477 


0.96 


NGC6284 


H 


3.13 ±0.04 


2.33 ± 0.04 


1.31 ±0.03 


-1.32 


0.87 


3.9 


- 


4.279 


0.97 


NGC6304 


L 


2.41 ±0.05 


1.47 ±0.04 


2.63 ± 0.03 


-0.59 


- 


0.0 


-1.00 


3.724 


1.51 


NGC6316 


L 


2.47 ±0.04 


1.46 ± 0.04 


2.13 ±0.02 


-0.55 


- 


- 


-1.00 


- 


- 


NGC6342 


H 


2.92 ± 0.11 


2.07 ±0.09 


2.00 ± 0.06 


-0.65 


0.94 


- 


-1.00 


3.778 


2.16 


NGC6352 


L 


1.98 ±0.05 


1.33 ±0.04 


2.37 ± 0.03 


-0.70 


- 


0.0 


-1.00 


- 


- 


NGC6356 


L 


2.40 ± 0.04 


1.58 ±0.03 


2.42 ± 0.02 


-0.50 


0.97 


0.0 


-1.00 


3.756 


1.12 


NGC6362 


L 


2.53 ±0.06 


1.99 ±0.05 


1.56 ±0.03 


-0.95 


0.92 


55.1 


-0.58 


3.954 


1.47 


NGC6388 


H 


2.89 ±0.02 


1.86 ± 0.02 


2.14 ±0.02 


-0.60 




2.4 


-0.70 s 


4.255 


1.01 


NGC6441 


H 


2.91 ±0.03 


1.80 ±0.02 


2.38 ± 0.02 


-0.53 




0.3 


-0.70 s 


4.230 


1.10 


NGC6528 


H 


3.00 ± 0.04 


1.59 ±0.03 


3.32 ± 0.02 


-0.04 




0.0 


-1.00 






NGC6553 


H 


2.94 ± 0.05 


1.52 ±0.04 


3.25 ± 0.03 


-0.21 




1.6 


-1.00 






NGC6569 


L 


2.53 ±0.08 


1.73 ± 0.07 


1.79 ±0.04 


-0.86 




0.0 




3.954 


1.12 


NGC6624 


L 


2.55 ±0.03 


1.66 ± 0.02 


2.23 ± 0.02 


-0.44 


0.88 


0.0 


-1.00 


3.771 


1.73 


NGC6637 


L 


2.35 ±0.03 


1.55 ±0.03 


2.11 ±0.02 


-0.70 


0.91 


0.0 


-1.00 


3.748 


1.41 


NGC6638 


L 


2.54 ± 0.05 


1.77 ±0.04 


1.79 ±0.03 


-0.99 


0.87 


18.3 


-0.30 


4.097 


1.14 


NGC6652 


H 


2.86 ± 0.03 


2.03 ± 0.02 


1.73 ±0.02 


-0.96 


0.92 


0.0 


-1.00 


4.000 


2.14 


NGC6723 


L 


2.64 ±0.05 


2.13 ±0.04 


1.26 ±0.03 


-1.12 


0.97 


20.5 


-0.08 


4.130 


1.05 



a GGCs with high (H) and low (L) H/3 indices at fixed metallicity. 

b Measured at FWHM = 3.lA spectral resolution. 

c From the Harris (1996, February 2003 revision) catalogue. 

d Relative GGC ages from De Angeli et al. (2005) and Recio-Blanco et al. (2006). 
° Maximum T c ff along the HB, from Recio-Blanco ct al. (2006). 

f ^BSS = ^°3(-^bss/-^F555w) inside the core radius. Lf555 w in units of 1O 4 L0. Taken from Moretti et al. (2008). 
s NGC6388 value taken from Zoccali et al. (2000). The same value is assumed for NGC6441 (Puzia et al. 2002). 
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Fig. 1. — H/3 (a) and H/3lick (b) indices versus [MgFe] for the 41 GGCs in S05. Solid and open symbols indicate, respectively, GGCs 
with high (GGCif) and low (GGCL) H/3 indices at fixed [MgFe]. SSP models from V99+ at FWHM = 3.lA are overplottcd, with dotted 
and dashed lines indicating fixed ages and metallicities as in the labels. 
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Fig. 2. — Metallicity versus CMD-derived parameters (AgeiMORMi SRRLyr) HBR, log(T a g HB), and Sggg; see definition in the text) for a 
subsample of 25 GGCs from S05 with [Fc/H]> —1.35. Symbol codes are kept as in Fig[T]except for the second parameter GGCs NGC6388 
and NGC6441 (open stars). Solid lines in panel e illustrate individual linear fits to solid and open squares. 



